Applications of the Wulff construction to the 

number theory 

Senya Shlosman* 
Centre de Physique Theorique, 
CNRS, Luminy, case 907, F-13288, 
Marseille Cedex 9, France, and 
IITP, RAS, Moscow, Russia 
shlosman@cpt.univ-mrs.fr 

February 7, 2008 



Abstract 

We apply the geometric construction of solutions of some varia- 
tional problems of combinatorics to estimate the number of partitions 
and of plane partitions of an integer. 
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1 Statement of results. 

In this note we will explain how certain classical statements from the number 
theory can be rederived using recent results of the asymptotic combinatorics. 
Specifically, we will derive the asymptotic behavior of the number of the 
partitions and of the plain partitions of an integer A^. 

We recall that a partition p of an integer is an array of non-negative 
integers rii > n2 > ... > > such that Xli^i — ^- One corresponds to 
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p a well-known geometric object called Young diagram, Y (p) , with cells. 
The set of all Young diagrams with cells will be denoted by y^. 

Similarly, a plane partition S of an integer is a two-dimensional array of 
non- negative integers Uij, such that for any i we have nn > ni2 > ... > riik > 

for any j we have nij > n2j > ... > rikj > while again YlTj=i ^ij ~ ^■ 
The corresponding geometric picture is called a skyscraper, or a mausoleum, 
or a 3D Young diagram. The set of all skyscrapers with cells (= the set 
of all plane partitions of A^) will be denoted by Sn- 

Our goal is to explain that as A^ — ^ oo 

2x1/2 



ln|3;;v|~7r(-) N'/\ (1) 



ln\Sj,\r^3{^Y' N'/'. (2) 



and 



In fact, much more precise information about the behavior of these functions 
is known. For |3^iv| it is the famous result of Hardy- Ramanuj an- Rademacher, 
An||. The corresponding asymptotic behavior of |iSAr| was found by 



see 



Wright in |[Wr| 



To get the relations (||) and @ , we are using, correspondingly, the results 
of ||VKer|| and [ pKe|| , where the shape of the typical Young diagram, resp. 
skyscraper, were found. (Here "typical" means with respect to the uniform 
distribution on the sets and S^.) These shapes are solutions of certain 
variational maximizing problems, described in details in the next section. 
Briefly speaking, these problems consist in finding a hypersurface in a 
certain class, which maximizes the surface integral (G) = j]jr]{nx) dsx, 
where n^^ is the normal to G at x, and is a corresponding given function. 

For the case of Young diagrams the function rjy is defined on a unit 
quarter-circle = S*^ fl M^, and is given by 

f \ 1 "1 , 1 
r]y (n) = — I rii m h n2 In ■ 



rii + n2 ni + n2 

The typical Young diagram with A^ cells, scaled by the factor A^~^/^, goes 
along the curve C : 

exp < —X > + exp < —y 
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This result was obtained by Vershik, see | VKer ] or [|V|. 

In the skyscraper case the function rj^ is defined on the set = S*^ fl M.^ 
of unit vectors n and is given by 



Vs 



1=1 \ I 11/ 



where |n|-^ = ni + 7^2 + t?-3, and L is the Lobachevsky function, defined for 
X G [0, tt] by 



Lix) 



In (2sint) dt. 



The surface in , which describes the shape of a typical skyscraper with 
cells scaled by the factor A^~^/^ to the unit volume, is given by 



G 



vs 



C(3) 



-1/3 



{ (/ {A, B,C) - In A, f {A, B,C) - In BJ {A, B, C) - In C) } 



where A, S, C > 0, A + B + C = 1, and 



f{A,B,C) = ^[ [ In {A + Be'"" + Ce''') dudv, 

47r J[o,27r] J[0,27t] 



see 



The geometric construction, which gives the solution to problems de- 
scribed above, was found in [^. In fact, it is very close to the Wulff con- 
struction, 1^]. Both constructions are described below, together with the 
relation between them. Due to the relation (|TT1), the evaluation of the inte- 
grals Vn (Grj) becomes very easy. It turns out that 



Vr^y (Gr^y) = TV ( - j , 



(3) 



while 



V,, {G 



'ns/ 



C(3) 



1/3 



(4) 



3 



We know from [ VKer | that most of the Young diagrams with cells go along 
the curve N^/'^C. On the other hand, if G is a graph of a decaying function 
g (x) with g (x) dx = N, then the meaning of the functional Vriy (G) is the 
following: "the number of the diagrams with cells going along the curve G 
is of the order of exp { V,;^ (G*) } " • Therefore, the value (||), multiplied by iV^/^ 
is precisely the exponent of (|l|). Correspondingly, most of the skyscrapers 
with cells go along the surface N^^^Gn^, 



CKe|. 



V,,{N'/'G,,)=N'/-'V,AGvs: 



the reader should consult the papers |[VKei 
paper 



see 

has to be the exponent 



Therefore the value 
1^). For more details 
|CKe|| and the review 



2 The maximizing problem. 

Let S"^ denote the d-dimensional unit sphere in M'^"'"^. Introduce the subset 
_ f-^-^drj-i ^£ "positive" unit vectors, lying in the positive octant. Let a 
function 77 on be given, which is supposed to be continuous, nonnegative: 
?7 (■) > 0, and vanishing on dA'^. We suppose for simplicity that the function 
1] vanishes also in some neighborhood of dA'^. (The general result which is 
needed for the formulas in the previous section is then obtained by letting 
this neighborhood to shrink to zero.) Let now G C M"^"*"^ be an embedded 
hypersurface, possibly with a boundary. We assume that for almost every 
xgG the normal vector ria, to G is defined, and moreover 

G A'^ for a.e. xeG. (5) 

Then we can define the functional 

Jg 

here ds is the usual volume (i-form induced from the Riemannian metric on 
M'^+i by the embedding G C R'^^^ 

Denote by Qat C M."^^ the cube, consisting of points x G R*^"*"^ with 
< Xi < N, i = l,2,...,d+l. Denote by O the point (0, 0, 0) G Qn, and 
let A = {N, N, ...,N) G Qn be the opposite vertex. 

Let the number V E (O, N'^'^^) be given. We introduce the family Dy of 
hypersurfaces G C Qat as follows: 

i) G splits the cube Qn into two parts, with the points O and A belonging 
to different parts; denote them hj Qn (G*, O) and Qn {G, A) ; 
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a) the property holds 

in) for the volume vol (G, A^) , defined as the {{d + l)-dimensional) volume 
vol {Qn {G, O)) of the body Qn (G, O) , we have 

vol (G, N) = V. 

We want to solve the following variational problem: find the upper bound 
of Vri over Dy : 

Vr, = sup (G) , (6) 

as well as the maximizing surface(s) ^ -Dy , such that V,, (K?) = ''^t?, if the 
maximizer does exist. 

It turns out that there exists a geometric construction, which provides a 
solution to the variational problem It was found in [ p^ . 

Let 

K> = {x G : Vn G A'^ (x, n) > r] (n)} , (7) 

G, = a (ir>) . (8) 

Note that is a convex set, so the hypersurface G^ has normals at almost 
every point. Moreover, these normals fall into A'^. Define the dilatation 
parameter A {V, N) as a unique solution of the equation 

vol {Qn n (Mf 1 \ A (y, N) K>)) = V. (9) 

Theorem 1 The unique solution to the variational problem ^ is the surface 

Gr,,N,v = QNn\{V,N)G^, (10) 
provided N is large enough. 

Remark 2 The analog of the variational problem which is obtained by 
the removal of the constraint is ill posed. 
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Note that the set is a convex region in M'^"'"^. If at some boundary 
point the region has a unique support plane, then this plane is of the 
form 

L^(n) = {xGM''+^:(x,n) = r/(n)} 
for the corresponding n. Therefore for large we have 

™,(G,)^1|(^. (11) 

The proof of the above theorem relies on the following construction, which 
corresponds to the variational problem (|^) a certain Wulff variational prob- 
lem. Once the correspondence is established, the proof follows easily from 
the known properties of the Wulff problem. 



3 The corresponding WulfF problem. 

First we formulate the general Wulff problem. We restrict ourselves to the 
symmetric case, in order to make the relation with the above maximizing 
problem more transparent. Let the real function r on S*^ be given. We 
suppose that the function is continuous, positive: r (■) > const > 0, and 
symmetric with respect to the reflections in coordinate planes: r (n) = 
r (ni, n2, rid+i) = r (±ni, ±n2, in^+i) for any choice of signs. Then 
for every hypersurface M'^ (possibly with a boundary), embedded in M'^"'"^, 
we can define the Wulff functional 

Wr{M^) = [ T(n,) ds,. (12) 

Here x G M'^ is a point on the manifold M'^, and the vector n^, is the unit 
vector parallel to the normal to M'^ at x. We suppose that Af^ is smooth 
enough. 

Suppose additionally that the hypersurface M'^ lies in fact in M.'^^, and 
separates the origin O from infinity in R++^ Let iV'^+^ C M++^ be the union 
of all finite components of W^'^ \ M'^; we denote the volume |A^'^"'"^| of A^'^^^ 
by vol (M*^) , and will call it the volume inside iVf^. We denote by D' the 
collection of all such hypersurfaces M*^ in with finite volumes. By 

D'^ C D' we denote the collection of all these hypersurfaces M'^, for which 
the volume vol (M'') inside M'^ equals A. 
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The Wulff problem consists in finding the lower bound of Wr over : 

Wr= inf W,(M), (13) 

as well as the minimizing surface(s) Wt^\ such that Wr {Wr^^^ = w,-, if it 
exists. 

The answer is given by the following Wulff construction. Let 

K< = {x G W^^^ : Vn (x, n) < r (n)} , and = M++^ H dK<. (14) 

Note that the set is a symmetric convex bounded region in M'^^^. If at 
some boundary point the region has a unique support plane, then this 
plane is of the form 

L,(n) = {xeM^+^ : (x,n) =T(n)} (15) 
for the corresponding n. Therefore 

VOl{Wr) = -^Wr{Wr). (16) 

d + 1 



Theorem 3 The variational problem ( [7^ has a unique solution 



Wr, 



which is a scaled version of the Wulff shape Wr, see (\T4). 



Remark 4 The standard use of the name "Wulff shape" refers to the full 
surface dK^ itself. In our symmetric setting the surface dK^ is a union of 
Wr and its multiple reflections in coordinate planes. 



Remark 5 The variational problem ( |7^ might have no solutions if the func- 
tion T is allowed to vanish. 
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The paper | [1'2|| contains a simple proof that Wr (^t^^^ < Wr (M) for 
every M G D'^. The uniqueness of the minimizing surface is proven in ||T1|| . 



Now we are going to construct the Wulff problem, corresponding to the 
problem (H) (and its solution (|T0|) ). For this we have to specify the function 
(n) and the value A of the volume. First we put 

A = N'^-V. 

We define r (n) under assumption that the normalization constant A {V, N) 
from the equation equals to 1. This is done to simplify the notations. In 
this case for n G A'^ we put 

(n) = dist (A, (n)) , 
see ([T5|). For the remaining n-s the function r is defined by symmetry. 
Note that 

Tr, (n) +r]{n)= (n, Ol^ . 
Therefore for G in Dy we have: 



V, (G) + Wr, (G) = (n,., Ol^ ds,. 



Let n be the hyperplane orthogonal to the vector oA, and P denote the 
orthogonal projection on 11. Then the last integral is nothing else but the 
(i-dimensional area S (P (G)) of the set P (G) , so 

Vr,{G)+Wr,{G)=S{P{G)), (17) 

and therefore 



sup Vrj{G)< sup S{PiG))- inf Wr, (G) (18) 

GeD^ G&D^ GeD^ 

Denote by Gq the surface made from d faces of the cube Qjy containing the 
vertex O. Then clearly S {P (Gq)) > S {P (G)) for every G G Dy . Since we 
suppose that the function r] vanishes in some neighborhood of c^A'^, we have 
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that S {P (Gfj)) = S {P (Go)) ) provided is large enough. Also, by the very 
definition of the function and in view of the Theorem ^ 

inf Wr, (G) = Wr, (G,) , 

SO ( p!7|) and (|18D imply that indeed 

sup Vrj (G) = (G^) . 
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